The importance of precession in modelling the 
direction of the final spin from a black-hole merger 



Enrico Barausse 

Center for Fundamental Physics, University of Maryland, College Park, MD 20742-4111, USA 

Abstract. The prediction of the spin of the black hole resulting from the merger of a 
generic black-hole binary system is of great importance to study the cosmological evolution 
of supermassive black holes. Several attempts have been recently made to model the spin 
via simple expressions exploiting the results of numerical-relativity simulations. Here, I first 
review the derivation of a formula, proposed in Ref. [T], which accurately predicts the final spin 
magnitude and direction when applied to binaries with separations of hundred or thousands of 
gravitational radii. This makes my formula particularly suitable for cosmological merger-trees 
and N-body simulations, which provide the spins and angular momentum of the two black holes 
when their separation is of thousands of gravitational radii. More importantly, I investigate 
the physical reason behind the good agreement between my formula and numerical relativity 
simulations, and nail it down to the fact that my formula takes into account the post-Newtonian 
precession of the spins and angular momentum in a consistent manner. 



The dynamics of black-hole (BH) binaries is a very complex problem which has been solved 
only very recently through time-expensive numerical-relativity (NR) calculations. However, 
it has been shown that the dimensionless spin of the remnant from a BH binary merger, 
ag n = Sfi n /M| n , can be described with simple prescriptions based on point particles [31 IU [5], 
on fits to the NR data [6J [9j QH] , or on a combination of the two approaches [11] . These 
formulas are useful because they provide information over the entire 7-dimensional space of 
parameters for BH binaries in quasi-circular orbits, namely: the mass ratio q = MijM\ and 
the six components of the initial dimensionless spin vectors a\ t 2 = S\ 2/^12- Such parameter 
space could in principle be investigated entirely via NR calculations; in practice, however, the 
simulations are very expensive and restricted to q = 0.1-1. Also, these formulas have applications 
in astrophysics, where they could provide information on the properties of massive-star binary 
systems; in cosmology, where supermassive BHs (SMBHs) are believed to assemble through 
accretion and mergers; in gravitational- wave astronomy, where the a priori knowledge of the 
final spin can help the detection. 

While the different expressions for the spin norm, |afl n |, are in good agreement with the 
results of NR simulations, the predictions for the final spin direction, ag n = afi n /|ctnn|) do 
not agree well with one another and are all essentially imprecise when the binaries are widely 
separated. This is because all expressions are built from and model the typical NR binaries 
and hence the dynamics of the last few orbits before the merger. Because it does not account 
systematically for the precession of the orbital angular momentum L, the prediction for a^ n 
depends on the separation of the binary and is therefore of little use for applications, such 
as cosmological merger-trees or N-body simulations, that provide the spins of the two BHs at 
separations of thousands of gravitational radii. Although one could use the PN equations to 



evolve a widely-separated binary to a separation of few gravitational radii and then apply the 
formulas, this makes the formulas difficult to use and implement. In Ref. pQ we followed instead 
a different approach and showed that it is possible to derive a formula for the final spin that 
takes into account the precession of the spins and that is, therefore, applicable to binaries with 
arbitrary separations. 

The derivation of the formula, for more details about which we refer the reader to Ref. pQ, is 
based on the following assumptions: 

(i) When the spins are parallel to L (either aligned or antialigned) , the formula must reduce to 
the fit of the NR results presented in Refs. fJJ, [77] /. namely 

afin = a + av(s±a + s$v + t ) + z/(2y3 + t 2 v + t 3 v 2 ) , (1) 
where v = M±M 2 /(Mi + M 2 ) 2 is the symmetric mass ratio, a = (a\ + 02a 2 ) /(l + q 2 ), ancQ 
s A = -0.1229 ± 0.0075 , s 5 = 0.4537 ± 0.1463 , 

to = -2.8904 ± 0.0359 , t 2 = -3.5171 ± 0.1208 , t 3 = 2.5763 ± 0.4833 . (2) 

(ii) The mass M rac j radiated to gravitational waves can be neglected i.e., Mfi n = M = M\ + M 2 . 
The reason why assumption (i) is reasonable here is that M rac j is largest for aligned binaries but 
these are also the ones fitted by expression ([jrj. In this way, the mass losses are automatically 
accounted for by the values of the coefficients to, t 2 , £3, S4 and S5. 

(Hi) The norms \S\\, \S 2 \, \£\ do not depend on the separation of the binary r, with £ being 

£(r) = S &n - [S^r) + S 2 (r)] = L(r) - J rad (r) , (3) 

where S\(r), S 2 (r) and L(r) are the spins and the orbital angular momentum at the separation 
r and Jrad(0 is the angular momentum radiated from r to the merger. While the constancy of 
\Si\ and \S 2 \ is a very good assumption for BHs, the constancy of \£\ is heuristic and based on 
the idea that the merger takes place at an "effective" innermost stable circular orbit (ISCO), so 
that \£\ can be interpreted as the residual orbital angular momentum contributing to S& n . 

(iv) The final spin S^ n is parallel to the initial total angular momentum J{r- m ) = S'i(rin) + 
S 2 (n n ) + L(r in ). 

(v) The angle between L and S = S\ + S 2 and the angle between S\ and S± are constant during 
the inspiral, although L and S precess around J. 

Assumptions (iv) and (v) are motivated by the PN approximation. It has been in fact 
shown by Ref. [13] that within the adiabatic approximation the secular angular-momentum losses 
via gravitational radiation are along J. This is because as L rotates around J, the emission 
orthogonal to J averages out. Therefore, S and L essentially precess around the direction J, 
which remains nearly constant (cf. the detailed discussion in Ref. j!3|). and the angle between 
L and S and that between the two spins remain constant as well. 

(vi) When the initial spin vectors are equal and opposite and the masses are equal, the spin of 
the final BH is the same as for nonspinning binaries. Besides being physically reasonable - if the 
spins are equal and opposite, their contributions are expected to cancel out - this assumption is 
confirmed by NR simulations and by the leading-order PN spin-spin and spin-orbit couplings. 

1 The number of free parameters of the fit is actually four, since to, £2, £3, S4 and S5 must satisfy the constraint 

afin = ^ + ^ + ll = 0.68646 ± 0.00004 . 
2 16 64 



which follows from the results obtained by Ref. [12] for equal-mass non-spinning BHs. 



These assumptions are sufficient to derive an expression for both the magnitude and the 
direction of the final spin pQ. In particular, the final spin norm is given by 



l a fm| - jp^yl 



a-i\ 2 + 1*3.2 1 2 ^ 4 + 2|a 2 ||ai|(7 2 cosa + 



where 
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2(|a 1 |cos/3+ |a 2 |g 2 cos 7 ) \£\q + \£\ 2 q 2 , (4) 
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and the angles a, (3 and 7 are defined by 

cos a = ai(r in ) • a 2 (r in ) , cos/3 = ai(r in ) • L(r in ) , cos 7 = a 2 (r in ) • L(r in ) . (5) 
The angle 9a n between the final spin and the initial orbital angular momentum L(r- m ) is instead 

cos 6> fin = L(r in ) ■ J (r in ) . (6) 

We will now test our expressions and ([5]) for | a^ Q | and our expression ([6]) for 6^ n against the 
NR simulations for generic binaries (i.e., with spins not parallel to L) published so far. Also, 
we will compare our predictions (AEI) with those of similar formulas suggested by Refs. [I] 
(BKL), [TT] (AEI old) and [8] (FAU). The comparison consists of two steps. First, we compare 
the different predictions using as input the initial data of the NR simulations, in which the 
binaries have small separations [r- in < 10 M). Second, using binaries at large separations 
( r in < 2 x 10 4 M), for which the dynamics starts being dominated by gravitational- wave emission 
and thus of direct relevance for cosmological investigations. More precisely, we evolve the NR 
initial configurations back in time up to a separation of 2 x 10 4 M using the PN equations in 
the quasi-circular limitd, calculating the predictions of the different formulas at each step and 
considering the maximum error for each formula. (The results that we present, however, do not 
change significantly if we integrate only up to a separation of ~ 200M.) 

The upper left panel of Fig. [1] shows the predictions of the various formulas for |ofi n |, 
when applied to the small-separation configurations corresponding to the initial data of the 
NR simulations (see caption for details). In particular, it reports the error ||afi n) NR,| — | a fm*||) 
where "*" stands either for "AEI" (which, as mentioned in Ref. [1], gives the same predictions 
for I cifin I as "AEI old"), "FAU" or "BKL". The lower left panel shows instead the maximum 
error when the configurations are evolved back in time up to r- m = 2 x 10 4 M. Although the AEI 
expression is slightly better, all the formulas give accurate predictions for the final spin norm, 
both for small and large separations. Note that the larger errors for the indices 1-10, which 
correspond to the simulations of Ref. [2] with small mass ratios (q = 0.13-0.17), are most likely 
due to the larger truncation errors affecting those simulations (see Ref. [13], sec. IIIA). 

The situation is very different when considering the final spin direction. In particular, the 
right panels of Fig. Q] report the inclination-angle error, |#fi nj * — arccos [L(ri n ) ■ afi ni NR,]|, for the 
data in the left panels, but for those of Ref. [T7], for which the final spin direction was not 
published. Clearly, when considering small-separation binaries (upper right panel), our new 
formula performs slightly better than the "BKL" and "AEI old" formulas, but it is not better 

2 In particular, following Ref. [2], we use the precession equations for the spins and the angular momentum at 2 
PN order and the rate of change of the frequency at 3.5 PN order (with spin terms included up to 2 PN order). 
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Figure 1. Left: The upper panel shows the error ||afi ni NR,| — |afi n .*|| ("*" being either "AEI" , 
"FAU" or "BKL" ) of the various formulas for the final spin norm, when applied to the small-separation 
configurations corresponding to the initial data of the simulations of Refs. [Tl] (indices 1-40), [TS] (indices 
41-42), [T5] (index 43), [5] (indices 44-76) and [T7] (indices 77-84). The lower panel shows instead the 
maximum error ||afi n ,NR| — |ofm.*| | when the configurations are evolved back in time up to large separations 
of r = 2 x 10 4 Af . Although the new AEI expression is slightly better, all formulas give accurate predictions 
for |ag n |, both for small and large separations. Note that the larger errors for the indices 1-10 are due 
to the larger truncation errors affecting those simulations (see text for details). Right: The same as in 
the left panel but for the inclination angle error |0fi n ,* — arccos [L(r in ) ■ afi nj NR]| and without the data of 
Ref. [T7], for which the final spin direction has not been published. The new AEI expression is accurate 
both for small and large separations, while the other ones become imprecise for large separations. 



than the "FAU" formula. Indeed, the latter is exact by construction for the indices 34-66. 
This is because for such data the final spin direction has not been published and it has been 
here reconstructed using the FAU formula applied to the configurations of Table II in Ref. [8]. 
However, when considering large-separation binaries (lower right panel), our new formula clearly 
performs much better than all the other ones. For instance, our error for 0^ n is below 7 degrees, 
for any separation, while it can be as large as 70 degrees with the older formulas. (The "steps" 
in the lower right panel reflect the different sequences of Table I of Ref. |14j.) 

To highlight the role played by the precession of the spins in the prediction of the final spin 
direction, in the left panel of Fig. [2] we focus on the binary "SP3" of Ref. [15] (index 41 in Fig. [I]). 
In this configuration the spins precess strongly and the final spin flips relative to S{r- m ). In 
particular, S\ and S2 are initially parallel, orthogonal to L and have: a\ = ai ~ 0.5, q = 1. The 
filled hexagon in Fig. [2]is the numerical result as obtained with an initial separation r- m ~ 6.6 M, 
while different lines show the angle #fi n obtained with the various formulas as a function of the 
binary separation. More specifically, using the initial conditions of the NR simulation we have 
evolved the "SP3" binary back in time using the PN equations up to a separation of 2 x 10 4 M, 
applying the various formulas set each, step to compute $fm,* 

= arccos [L(r- m ) ■ 0^*]. The line 
"PN+NR" reports instead the angle arccos [£<(ri n ) • ag^NR,]- In the right panel, we consider the 
configuration "Q13TH00" of Ref . [2] (index 1 in Fig. [TJ. This configuration features a smaller, 
non-spinning BH and a larger, spinning BH with spin a ~ 0.81 lying on the equatorial plane and 
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Figure 2. Left: Predictions of the different formulas for 9f\ n for the binary "SP3" of Ref. [TS] as a 
function of the initial separation r; n . Shown with the "NR+PN" curve is the angle arccos [L(r- m ) ■ S^nr], 
where L(r- m ) is obtained by integrating the PN equations starting from the initial conditions of the NR 
simulation, while the filled hexagon shows #fi n as computed by the NR simulation. Right: The same as 
in the left panel, but for the configuration "Q13TH00" of Ref. [33]. Clearly, only our formula reproduces 
the "NR+PN" results to within 2-3 degrees for any r ln , while the other formulas have errors that grow 
very large for separations r m > 200M. 



orthogonal to the separation (r; n ~ QAM). The mass ratio is q ~ 0.13, which makes this binary 
interesting as it allows to test the behavior of the various formulas in a region of parameter 
space which so far has been sampled only sparsely by NR simulations, which usually deal with 
q ~ 1 . (This is because the cost of a NR simulation grows as 1/q 2 , a factor 1 jq coming from the 
slow inspiral of a small mass ratio binary, and another factor 1/q appearing because if the mass 
ratio is small, one needs much resolution to resolve the two very different scales of the problem.) 

Three results are clear from these comparisons. First: all previous formulas provide a 
reasonable estimate of 9a n but only when the binary has a separation r- m <C 200Af . Second: 
Only our formula provides a reasonable estimate at all separations; indeed, the "PN+NR" 
results are reproduced to within the accuracy of the NR simulation, i.e. 2-3 degrees. All the 
other formulas, instead, predict the same value for 6^ n for all r; n , because they take as an input 
the angles arccos(5i ) 2 • L) and arccos(S'i • S2), which are constant under the quasi-circular PN 
evolution |2j. Therefore, the other formulas become rapidly become imprecise when the binary's 
separation increases, reaching errors as large as 60 degrees for r- m > 200M . Third: the error of 
the old formulas when applied to large separation binaries is larger for small mass ratio systems 
{e.g. the maximum error is ~ 13 degree for the comparable-mass configuration SP3, while it is 
~ 60 degrees for the small mass-ratio binary Q13TH00). This had to be expected, again on 
the grounds that the old formulas predict the same #fi n for all r; n . As a result of this, in fact, 
the maximum error of the old formulas is roughly given by their prediction for 8^ Q at small 
separations, because the correct 9^ n becomes small for large r- m (c/. the "PN+NR" curves in 
Fig- G])- This prediction can be very large for small q if the angle between S and L at small 
separations is large. This can be easily understood by noting that at small separations the old 
formulas give roughly the same results as our formula: Therefore, for q ~ their prediction is 



cos(0 fin ) « L(r- m ) ■ S(r in ) [because for q 0, J(r in ) « S^in)]- 

We stress that such a good agreement between our formula and the data for the final spin, 
irrespective of the separation of the binary to which the formula is applied, emerges because 
we have consistently taken into account the effect of precession through assumptions (iv) and 
(v), while that effect is not accounted for in the other formulas. For instance, in our earlier 
formula of Ref. [11] ( "AEIold" ) , assumptions (iv) and (v) were replaced by the assumption that 
the angular momentum emitted during the inspiral is parallel to L(r- in ) [cf. assumption (Hi) 
of Ref. [H]]. This is not true unless one neglects spin-orbit precession (as was indeed stressed 
explicitly in Ref. [TT] ) . A similar assumption was also made in Ref. [I] , while Ref. [T7] admittedly 
recognized that their formula might not be valid at arbitrary separations, as it is based on a fit 
to NR simulations, which, as already mentioned, have initial separations r- m ~ 10M. 

In conclusion: I have reviewed the assumptions needed to derive a new formula predicting 
the spin of the BH resulting from the merger of two BHs in quasi-circular orbits and having 
arbitrary initial masses and spins [I]. This formula includes the effect of the precession of 
the spins through assumptions (iv) and (v), and can therefore be applied to widely separated 
binaries, such as those relevant for cosmological applications, for which the other available 
formulas become imprecise. I stress that requiring that the formulas for the final spin work 
also when applied to large separation binaries is necessary to ensure that these formulas can 
be readily used in cosmological applications. Cosmological simulations typically provide the 
initial spins and angular momentum of SMBH binaries when their separation is still of about 
0.1 pc (see Ref. [18] for a recent example), which corresponds to 2 x 10 5 Af if M = 10 8 Mq. 
This separation roughly corresponds to the point at which the evolution of the binary starts 
being driven by gravitational-wave emission. While it would be in principle possible to use the 
PN equations and evolve the widely-separated binaries provided by cosmological simulations to 
small separations, read-off the relevant information and apply the formulas for the final spin, 
this procedure would be impractical and potentially very time-consuming (cf. again Ref. |18j . 
which simulates tens of thousands of SMBH binaries in order to get significant statistics). 
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